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Abstract 

Recently the DAMA/Nal Collaboration reported further indication of a 
possible modulation effect in WIMP direct detection. In this note we discuss 
the relevance of this result for supersymmetric theories where the neutralino 
is the dark matter candidate. We specifically consider the Minimal Supersym- 
metric extension of the Standard Model (MSSM) and Supergravity-inspired 
schemes, with possible deviations of the unification conditions at the GUT 
scale in the Higgs sector. The main results of our analysis show that the an- 
nual modulation data are widely compatible with an interpretation in terms 
of a relic neutralino as the major component of dark matter in the Universe. 
We also discuss the implications for searches of supersymmetry at accelerators 
and for indirect searches of neutralino dark matter. 
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Recently the DAMA/Nal Collaboration reported further indication of a possible 
modulation effect in WIMP direct detection. In this note we discuss the relevance 
of this result for supersymmctric theories where the neutralino is the dark matter 
candidate. We specifically consider the Minimal Supersymmctric extension of the 
Standard Model (MSSM) and Supcrgravity-inspircd schemes, with possible devia- 
tions of the unification conditions at the GUT scale in the Higgs sector. The main 
results of our analysis show that the annual modulation data are widely compat- 
ible with an interpretation in terms of a relic neutralino as the major component 
of dark matter in the Universe. We also discuss the implications for searches of 
supcrsymmetry at accelerators and for indirect searches of neutralino dark matter. 

1 Introduction 

The presence of non-baryonic dark matter in our Galaxy can be probed by 
means of different techniques which attempt to detect either directly (through 
elastic scattering off nuclei) or indirectly (through products of annihilation) 
the dark matter particles which are supposed to be embedded in the galactic 
halo. It has been shown (see, e.g., Ref. : ) that the sensitivities of the present 
experiments is currently at the level required for the study of one of the most 
appealing particle candidates for dark matter, the neutralino. This implies 
that it is now feasible to start the investigation of possible signatures in the 
detection rates, originated by specific features related to the presence of the 
dark matter particles. 

In the case of direct detection, a typical signature consists in the annual 
modulation of the detection rate. This effect was first pointed out in the semi- 
nal papers of Refs. 2 , where it was observed that, during the orbital motion of 
the Earth around the Sun, the change of direction of the relic particle veloci- 
ties with respect to the detector induces a time dependence in the differential 
detection rate, i.e. S(E, t) = Sq(E) + S m (E) cos[w(i — to)], where uj = 27r/365 
days and to = 153 days (June 2 nd ). So(E) is the average (unmodulated) dif- 
ferential rate and S m (E) is the modulation amplitude of the rate. The relative 
importance of S m (E) with respect to So(E) for a given detector, depends both 
on the mass of the dark matter particle and on the value of the recoil energy 
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where the effect is looked at. Typical values of S m {E) / Sq(E) for a Nal detec- 
tor range from a few percent up to ~ 15%, for WIMP masses of the order of 
20-80 GeV and recoil energies below 8-10 KeV. 

Over the last year, the DAMA/Nal Collaboration reported on two different 
analyses of the data collected during two periods of data taking, obtained with 
an experimental set-up consisting of nine 9.70 Kg Nal(Tl) detectors 3,4 . The 
data have been analysed by employing a maximum likelihood method, which 
allows to test the hypothesis of the presence of a yearly modulated signal 
against a time-independent background, by properly considering the energy 
and time behaviour expected for a recoil of massive WIMPs. The remarkable 
result of Refs. 3,4 is that the data support the possible presence of an annual 
modulation effect induced by WIMPS in the counting rate of the detector. By 
combining the data of the two periods of data taking, a total exposure of 19,511 
kg x day has been collected. The maximum likelihood analysis indicates that 
the hypothesis of presence of modulation against the hypothesis of absence of 
modulation is statistically favoured at 99.6% C.L., and pins down a 1-a C.L. 

region in the plane C^scaiar'™'' ~~ m x' wnere m x 1Sl ^ ne WIMP mass, f^".^ 001 ^ is 
the WIMP-nucleon scalar elastic cross section and £ = p x j pi is the fractional 
amount of local WIMP density p x with respect to the total local dark matter 
density pi. This region is plotted in Fig.l as a closed dashed curve, for pi = 0.3 
GeV cm~ 3 . The ensuing 1—a ranges for the two quantities are: m x = 59j^| 

GeV and fr^ilr"^ = 7.0±?; 4 x lO" 9 nb 4 . These results refer to v lms = 270 Km 
s _1 for the velocity dispersion of the WIMP Maxwellian velocity distribution 
in the halo, v csc — 650 Km s _1 for the WIMP escape velocity and u© = 232 
Km s _1 for the velocity of the Sun around the galactic centre. We notice that 
a variation of these velocity parameters inside their 1—a ranges of uncertainty 
affects the maximum likelihood value of m x and Cc/"^ 50 "- 1 by roughly ± 30% 
and by ~ ± 12%, respectively. 

In this paper, which is based on the results obtained in Refs. 5 ' 6 ' 7 we derive 
the theoretical implications of the experimental data of Refs. 3,4 , assuming that 
the indication of the possible annual modulation is due to relic ncutralinos and 
we show that these data are fully compatible with an interpretation in terms 
of a relic ncutralino as the major component of dark matter in the Universe. 
We select the relevant susy configurations in two different frameworks: the low 
energy minimal supersymmetric extension of the standard model (MSSM) and 
supergravity-inspired schemes (SUGRA models). In the latter case, we will 
explicitly consider both a strict unification at the GUT scale and the possibility 
of deviation from universality in the Higgs sector. The present analysis extends 
a previous analysis of ours 8 referring to the experimental data of Rcf. 3 . We 
then finally also discuss how the susy configurations, selected by the annual 
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Figure 1: Scatter plot of MSSM configurations compatible with the annual modulation data 

in the plane mx~€ a icala.r ■ The dashed contour line delimits the 2-cr C.L. region, obtained 
by the DAMA/Nal Collaboration, by combining together the data of the two running periods 
of the annual modulation experiment. The solid contour line is obtained from the dashed 
line, which refers to pi = 0.3 GeV cm~ 3 , by accounting for the uncertainty range of pi. (a) 
and (b) refer to configurations with 0.01 < f! x /i 2 < 0.7 and with U x h 2 < 0.01, respectively. 

modulation data, can be investigated by indirect searches for relic WIMPs and 
at accelerators. 

2 Selection of susy configurations by the annual modulation data 

In order to discuss which region in the susy parameter space is selected by 
the DAMA data, we first convert the region delimited by the 2-cr C.L. dashed 
contour line of Fig.l into an enlarged one, which accounts for the uncertainty 
in the value of pi , due to a possible flattening of the dark matter halo and a 
possibly sizeable baryonic contribution to the galactic dark matter: 0.1 GeV 
cm" 3 < pi < 0.7 GeV cm 3 . One then obtains the 2-cr C.L. region denoted by 
a solid contour in Fig.l (hereafter denoted as region R). The susy configura- 
tions arc then selected by the requirement that (m x , ^scalar*™' ') ^ R> when m x , 
^scalar an d £ are evaluated in the MSSM and SUGRA schemes. As for the 
values to be assigned to the quantity £ = p x /pi we adopt the standard rescaling 
recipe 5 : £ = min[l, f2 x /i 2 /(f2/i 2 ) m i n ], where Q, x h 2 denotes the neutralino relic 
abundance, calculated in the susy model 9 , and (flh 2 ) mm is a minimal value 
compatible with observational data and with large-scale structure calculations. 
We use here the value (f2/i 2 ) m i n = 0.01 5 . In all our analyses, we consider as 
cosmologically acceptable all configurations which provide £l x h 2 < 0.7. 
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Figure 2: Scatter plot of MSSM configurations compatible with the annual modulation 



data in the plane Q x h - £<j 



(nuclcon) 
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The two vertical solid lines delimit the Q x h -range of 



cosmological interest. The two dashed lines delimit the most appealing interval for Q x h , as 
suggested by the most recent observational data. The hatched area is excluded by cosmology. 



3 Analysis in the MSSM 

In our first analysis 5 we employ the minimal supersymmetric extension of the 
standard model (MSSM) which conveniently describes the susy phenomenology 
at the electroweak scale, without too strong theoretical assumptions. The 
MSSM is based on the same gauge group as the Standard Model, contains 
the susy extension of its particle content and two Higgs doublets. The free 
parameters of the model are: the SU(2) gaugino mass parameter M 2 , related to 
the U(l) gaugino mass parameter by the GUT relation M\ — (5/3) tan 2 9wM 2 ; 
the Higgs-mixing parameter /x; the ratio of the two Higgs vev's tan /3; the mass 
of the pseudoscalar neutral Higgs ttla; a common soft-mass parameter for all 
the squarks and sfermions mo; a common trilincar parameter for the third 
family A (the other trilinear parameters are all set to zero). The parameters 
are varied in the following ranges: 10 GeV < M 2 < 500 GeV, 10 GeV < \fj,\ < 
500 GeV, 75 GeV < m A < 1 TeV, 100 GeV < m < 1 TeV, -3 < A < 
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Figure 3: Explorability at accelerators of the MSSM configurations compatible with the 
annual modulation data, (a): Scatter plot in the plane - tan (i. The hatched region on 
the right is excluded by theory. The hatched region on the left is excluded by present LEP 
data at 183 GeV. The dotted and the dashed curves denote the reach of LEP2 at energies 
192 GeV and 200 GeV, respectively. The solid line represents the 95% C.L. bound reachable 
at LEP2, in case of non discovery of a neutral Higgs boson, (b): Scatter plot in the plane 
m x - tan/3. The hatched region on the left is excluded by present LEP data. The dashed 
and the solid vertical lines denote the reach of LEP2 and TeV33, respectively. 

+3, 1 < tan/3 < 50. 

Our susy parameter space is constrained by the latest data from LEP2 
on Higgs, neutralino, chargino and sfermion masses 11 . Moreover, the con- 
straints due to the b — > s + 7 process has been taken into account, considering 
the latest results both in the theoretical evaluation and in the experimental 
determination of the branching ratio 12 . 

By varying the susy parameters inside the ranges defined above, we find 
that a large portion of the modulation region R is indeed covered by susy 
configurations, compatible with all present physical constraints. This set of 
susy states (set S) is displayed in Fig.l with different symbols, depending 
on the neutralino composition. In Fig. 1(a) we notice that a quite sizeable 
portion of region R is populated by susy configurations with neutralino relic 
abundance inside the cosmologically interesting range 0.01 < fl x h 2 < 0.7. 
Thus we obtain the first main result of our analysis, i.e. the annual modulation 
region, singled out by the DAMA/Nal experiment, is largely compatible with 
a relic neutralino as the major component of dark matter. This is certainly 
the most remarkable possibility. However, we also keep under consideration 
neutralino configurations with a small contribution to Q x h 2 (see Fig. 1(b)), 
since also the detection of relic particles with these features would provide in 
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itself a very noticeable information. 

The neutralino relic abundance fl x h 2 is plotted versus S.o-^^ ^ in Fig. 2. 
We notice that a large fraction of the neutralino relic abundance falls into 
the restricted range 0.02 < ficDM^ 2 < 0.2, which turns out to be the most 
appealing interval for relic ncutralinos as indicated from recent observations 
and analyses on the value of the matter content of the Universe 5 . 

The properties of set S relevant to searches at accelerators arc displayed 
in Fig. 3. Section (a) of this figure shows a scatter plot of set S in term of 
m/j and tan/3, where it is apparent a correlation between tan/3 and m^. This 
is due to the fact that the rather large values Cg Ca i 

as required by the annual modulation data, impose that either the couplings 
are large (then large tan/3) and/or the process goes through the exchange of 
a light particle. Thus, Higgs-exchange dominance (which turns out to occur 
here) and ^",^ r con ^ ~ (10~ 9 — 10~ 8 ) nb require a very light h at small tan/3, 
and even put a lower bound on tan (3: tan/3 > 2.5. From Fig. 3(a) we notice 
that a good deal of susy configurations are explorable at LEP2, while others 
will require experimental investigation at a high luminosity Fcrmilab Tevatron, 
which should be capable to explore Higgs masses up to ~ 130 GcV. In 
Fig. 3(b) we display the scatter plot of set S in the plane m x - tan/3. Since the 
reach of LEP2 extends only up to the dashed vertical line, at m x ~ 50 GeV, 
the exploration of the whole interesting region will require Tevatron upgrades 
or LHC. Under favorable hypothesis, TeV33 could provide exploration up to 
the vertical solid line. 
4 Analysis in SUGRA schemes 

In this Section we show that the susy features, implied by the DAMA/Nal data, 
are also compatible with more ambitious supersymmetry schemes, where the 
previous phenomenological model is implemented in a supergravity (SUGRA) 
framework, especially if the unification conditions, which are frequently im- 
posed at the Grand Unification (GUT) scale, are appropriately relaxed 6 . 

The essential elements of the SUGRA models employed here are: a Yang- 
Mills Lagrangian, the superpotential and the soft-breaking Lagrangian. In this 
class of models the electroweak symmetry breaking (EWSB) is induced radia- 
tively. This supergravity framework is usually implemented by some restrictive 
assumptions about unification at Mqut'- (i) unification of the gaugino masses: 
Mi(Mcur) = m-i/21 (h) universality of the scalar masses: m,i{McuT) = mo, 
(iii) universality of the trilinear scalar couplings: A 1 {Mqut) = A d (McuT) = 
A m (Mout) = A m - As extensively discussed in Ref. 10 , these conditions 
have strong consequences for low-energy supersymmetry phenomenology, and 
in particular for the properties of the neutralino. 

The unification conditions represent a theoretically attractive possibility, 
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Figure 4: Scatter plot of SUGRA configurations compatible with the annual modulation 

data in the plane Q x h 2 - ^scalar*™^ . (a): universal SUGRA models, (b): SUGRA models 
with deviations from universality in the Higgs sector. Notations are as in Fig.2. 

which makes strictly universal SUGRA models very predictive. However, the 
above assumptions, particularly (ii) and (hi), are not fully justified, since uni- 
versality may occur at a scale higher than Mqut, he. the Planck scale or 
string scale, in which case renormalization above Mqut weakens universal- 
ity. We therefore discuss the DAMA/Nal data both in a SUGRA model with 
strict unification conditions and in a SUGRA framework, where we introduce 
a departure from universality in the scalar masses at Mqut which splits the 
Higgs mass parameters: M^XMqjjt) = ttlq(1 + Si). The parameters Si will 
be varied in the range (— 1,+1), but are taken to be independent of the other 
susy parameters. 

Because of the requirements of radiative EWSB and of the universality 
conditions, the independent susy parameters are reduced to (apart from the 
S^s): mi/2, m 0} A , tan/3 and sign(^). They are varied in the following ranges: 
10 GeV < mi/2 < 500 GcV, m a < 1 TeV, -3 < A a < +3, 1 < tan/3 < 50; 
the parameter \x is taken positive. The values taken as upper limits of the 
ranges for m 1 / 2 ,mo are inspired by the upper bounds which may be derived 
for these quantities in SUGRA theories, when one requires that the EWSB, 
radiatively induced by the soft supersymmetry breaking, does not occur with 
excessive fine tuning. The same argument was also used in the previous Section 
in setting the upper limits on the dimensional parameters of the MSSM. 

The susy parameter space is constrained by the same experimental bounds 
discussed in the previous Section for the MSSM, with the additional constraint 
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Figure 5: Explorability at accelerators of the non— universal SUGRA configurations com- 
patible with the annual modulation data, (a): Scatter plot in the plane mj, - tan/3, (b): 
Scatter plot in the plane m x - tan (i. Notations are as in Fig. 3. 



arising from the limits on the bottom-quark mass m&. The bottom mass is 
computed as a function of the susy parameters and required to be compatible 
with the present experimental bounds 13 . 

The susy configurations compatible with the DAMA data are shown in 
the plane ^cF^^ on ^ ~ ^ x ^ 2 m Fig. 4(a) for universal SUGRA models and in 
Fig. 4(b) for models with deviation from strict universality. We notice that also 
in SUGRA theories a fraction of the selected susy configurations fall into the 
cosmologically interesting range of Vl x h 2 . 

Other qualifications for the configurations which lie inside the region R, 
which are relevant for searches at accelerators, concern the ranges for the h— 
Higgs boson mass, the neutralino mass and the lightest top-squark mass. In 
the case of universal SUGRA models, we find: rrih < 115 GeV, 50 GeV < m x < 
100 GeV, 200 GeV < m~ u < 700 GeV and tan/3 > 42. For deviations from 
universality, the situation is shown in Fig. 5 where we notice that the sample of 
representative points covers a slightly wider domain. The ranges of the Higgs 
and neutralino masses are similar to those already found in the universal case, 
but now tan/3 extends to the interval 10 < tan/3 < 50, instead of being limited 
only to very large values. 

5 Indirect detection of neutralino dark matter 

The susy configurations which have been proved to be compatible with the 
annual modulation data can also be searched for by using methods of indirect 
search for relic particles . The two most promising techniques are the mea- 
surement of cosmic-ray antiprotons 14 and the measurement of neutrino fluxes 
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Figure 6: Scatter plots for the up-going muon fluxes from the center of the Earth versus the 
ncutralino mass. The MSSM configurations compatible with the annual modulation data 
are subdivided into the 4 panels, depending on the corresponding value the local density: 
p;/(GeVcm -3 ) = 0.1, 0.3, 0.5, 0.7. Dots denote configurations which could be excluded on 
the basis of the BESS95 antiproton data, circles denote configurations which survive this 
exclusion criterion. The dashed line denotes the MACRO upper bound. 

from Earth and Sun 15 . 

The signals we are going to discuss here consists of the fluxes of up- going 
muons in a neutrino telescope, generated by neutrinos which are produced by 
pair annihilations of neutralinos captured and accumulated inside the Earth 
and the Sun. The calculation of the up-going muon signal is performed ac- 
cording to the method described in Ref. 15 . 

In Fig. 6 we display the scatter plots for the flux of the up-going muons 
from the center of the Earth, for various values of the local dark matter den- 
sity p\. In this figure is also reported the current 90% C.L. experimental upper 
bound on $Earth^ obtained by MACRO 16 . We notice the particular enhance- 
ment at m x ~ (50 — 60) GeV, due to a mass-matching effect between the 
WIMP and the Fe nuclei in the Earth core. It is also noticeable an effect of 
suppression and spreading of the fluxes for m x > 70 GeV at pi > 0.3 GeV 
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cm~ 3 . This is due to the fact that the configurations with these large values 
of pi may imply, because of the annual modulation data, a ncutralino nucleus 
cross-section which is too small to establish an efficient capture rate, necessary 
for the capture- annihilation equilibrium in Earth 7 . In Fig. 6 the configura- 
tions which would be excluded on the basis of the antiproton data are denoted 
differently from those which would survive this criterion (for details see Refs. 

7,17)_ 

By comparing our scatter plots with the experimental MACRO upper 
limit, one notices that a number of susy configurations provide a flux in ex- 
cess of this experimental bound and might then be considered as excluded. 
However, it has to be recalled that a possible neutrino oscillation effect may 
be operative here and affect the indirect neutralino signal as well as the back- 
ground consisting in atmospheric neutrinos. Therefore a strict enforcement 
of the current upper bound on $Earth gjjQ^j applied with caution as long 
as the neutrino oscillation properties are not fully considered. However, it 
is rewarding that the set S of susy configurations is quite accessible to relic 
neutralino indirect search by measurements of up-going fluxes. 
6 Conclusions 

In this paper, we have analysed in terms of relic neutralinos the total sample 
of new and former DAMA/Nal data 3:4 , which provide the indication of a 
possible annual modulation effect in the rate for WIMP direct detection. The 
remarkable result of our analysis is that the annual modulation data are widely 
compatible with a relic neutralino making up the major part of dark matter in 
the Universe, both in the low energy MSSM and in SUGRA schemes. 

We have also investigated the possibility of exploring at accelerators the 
same susy configurations which are compatible with the annual modulation 
data. We have shown that an analysis of the main features of these susy 
configurations is within the reach of present or planned experimental set-ups. 
In particular, we have obtained the following results: 

• MSSM: The sizeable neutralino-nucleon elastic cross-sections, implied by 

the annual modulation data, entail a rather stringent upper bound for rah 
in terms of tan [3. In particular, this property implies that no susy config- 
uration would be allowed for tan/3 < 2.5. Another property, discussed in 
Ref. 5 , is that the annual modulation data and the b — > s + 7 constraint 
complement each other in providing a correlation between tan (3 and the 
mass of the lightest top-squark. 

• SUGRA: In the universal SUGRA model the constraints imposed by the 

DAMA/Nal data imply for the /i-Higgs boson mass, the neutralino mass 
and the lightest top-squark mass, the following ranges: mh < 115 GeV, 
50 GeV <m x < 100 GeV and 200 GeV < m fl < 700 GeV, respectively. 



10 



In universal SUGRA tan (3 is constrained to be large, tan/3 > 42, whereas, 
with departure from universality in the scalar masses, the range for tan (3 
widens to 10 < tan (3 < 50. 

Many of the above configurations will be explored by LEP2, and almost 
all of them are under reach of the future planned high-energy accelerators, 
namely the upgrade of the Tevatron and LHC. 

The same configurations can also be probed by indirect dark matter searches. 
We have shown 7 that a sizeable fraction of the susy neutralino configurations 
singled out by the DAMA/Nal data may provide signals detectable by mea- 
surement of cosmic-ray antiprotons and detection of neutrino fluxes from Earth 
and Sun. 

For the case of cosmic antiprotons, it has been shown 7 ' 17 that present data 
are well fitted by total spectra which include a p contribution from neutralino 
pair annihilation, with neutralino configurations which are relevant for annual 
modulation in direct detection. These data can also be used to reduce the total 
sample of the susy configurations under study, and to narrow the range of the 
local density, by disfavoring its largest values. Investigation by measurements 
of cosmic p's looks very promising in view of the collections and analyses of 
more statistically significant sets of data in the low-energy regime which are 
currently under way and which may soon provide further relevant information, 
like, for instance, in the case of BESS and AMS detectors. 

Measurements of neutrino fluxes from Earth and Sun, due to capture and 
annihilation of neutralinos inside these celestial bodies, have been proved to 
be sensitive to neutralino configurations singled out by the annual modulation 
data. However, an appropriate interpretation of these measurements prelim- 
inarily requires some clarification of the oscillation neutrino properties, espe- 
cially in the light of the recent Kamiokande result 18 . 

Finally, we conclude with a few comments. A solid confirmation of the 
annual modulation effect, as singled out by the DAMA/Nal Collaboration, 
necessarily requires further accumulation of data with very stable set-ups over 
a few years, a project which is currently undertaken by the DAMA/Nal Col- 
laboration. For the status of the experimental activity in WIMP direct search 
with other detectors, see the contributions of the various experimental Collab- 
orations to these Proceedings. 

It is really worth noticing that the detection of the effect of annual modu- 
lation, if confirmed by further experimental evidence, would turn out to be a 
major breakthrough in establishing the existence of particle dark matter in the 
Universe, and this would henceforth be a major breakthrough for astrophysics, 
cosmology and particle physics as well. 
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